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Abstract 

The fundamental solution (Green's function) of a first order matrix ordinary 
differential equation arising in a Landau-type problem is calculated by two methods. 
The coincidence of the two representations results in the integral formula for the 
product of two parabolic cylinder functions Dy{x)Dy{—x) at Re < 0, a; G iR. 
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1 INTRODUCTION 

This paper is to point out another integral representation for the product of two parabohc 
cyhnder functions D^{x)D^{—x) at Reu < 0,x E M. The formula to be obtained is due 
to equivalence of two representations for the fundamental solution (Green's function) of 
a first order (2 x 2-matrix) ordinary differential equation related to a couple of Weber's 
equations. The parabolic cylinder functions have been extensively studied in classical 
mathematical physics [1-5] (and refs. therein). In particular, among various relations 
for them, the following integral representations are known for products of two parabolic 
cylinder functions: 

D^z) D^,_i{z) , (1.1) 
{z real, Reu < 0, or Rez > O^Reu > —1), 

D_,_i(^e^^)L'_,_i(^e-^i), (1.2) 

(Rez > 0, Reu > -1, or \argz\ < f , Rei/ < -1, or -1 KReu < 0), 

D,{ze'^ D,{ze-'^^) , (1.3) 

{\argz\ < f,Rei/ < 0) [1, 2, 6]. Besides, the following relations can be found in [5]: 



and [6]: 



1 (7 f°° 

- (D,(z) D,{-z)) = Uzs)Dl{s)sds , (2) 
z az Jo 



{2z)-^D,{-z) D,+,{z) - D,+,{-z) D,{z)) 



{2zr' (z - D,{z)D,{-z) 



/■oo 

= / Jq{zs)D,{s) D,^^{s)ds . (3) 
Jo 

Moreover, Ref.[l] suggests to consult with [7, 8, 9] for other information concerning inte- 
gral representations for the parabolic cylinder functions. Essentially, apart from (1), only 
the following combinations can be found in [7-16]: 

(L'.(^e't) + D,{ze-'^^ (D_,_,{ze'^) + D_,_^{ze-^'^)) (4.1) 

(-1 < Reu < 0), 

{D,{ze'^2) - D,{ze-'l)) {d _,_,{ze'^^) - D_,_,{ze-'l)) (4.2) 
(-2 < Rei/ < 0), 

D,{z) {D,{ze'^) ± D,(ze-'^)) (4.3) 

{\argz\ < 7r/4, Rei/ < 0)). There are more other representations for products of D,y{z) by 
means of integrals and series in [6, 17-24]. For instance, the products at the same argu- 
ments can be found in [19, 21-23]. Certain representations through complex and contour 
integrals are given in [22]. Expressions for products by means of indefinite integrals are 
given in [19]. Recently, it has been found that the wave field of the plane wave, which is 



2 




scattered by a cone of arbitrary shape, can also be expressed for singular directions by 
means of the parabolic cyhnder functions (the scalar and electromagnetic cases) [25] . 

Combining (2) and (3) we obtain the following integral for the product of two functions 
Di,[z) at the opposite arguments: 

POO 

D,{z)D,{-z) = -2 / J^{zs)D,{s)Dl{s)ds (5) 

Jo 

(notice that integration of (2) agrees with (5)). In its turn, the present paper is to point 
out another formula for the product D^{x) D^{—x): 

D,{x)D,{-x) = 2-V2r-i(-z/) 

a;2 t\ dt 
- — tanh - , . , Rei'<0,xe]R, (6) 
2 ^ J V sinh t 

which appears as a by-product of [26, 27] where a theoretical problem of condensed mat- 
ter physics is studied. Specifically, the ground state average of momentum operator (so- 
called, mass current) in weakly inhomogeneous A- phase of superfluid helium-3 (^He— A) is 
investigated in [26, 27]. Equation (6) appears since the Green function of the quantum- 
statistical model under consideration (i.e., the fundamental solution to an appropriate 
matrix ordinary differential equation) can be obtained in two equivalent forms: the inte- 
gral [26], and the series [27]. Equating to each other, one obtains (6). The differential 
equation under consideration below reminds the famous physical problem of quantization 
for a spinning electron in a constant homogeneous magnetic field (the Landau problem) 
[28]. 

It is especially convenient to apply to (6) the Laplace method (steepest descent) in 
order to calculate asymptotical expansions for the mass current [27]. Therefore, Eq.(6) 
looks attractive since can be considered at the same footing as (l)-(5), and, hopefully, 
can be useful for applications, since contains only elementary functions under integration. 
The present paper briefly recalls the main points of [26] and [27] in Sections 2 and 3, 
respectively. Section 4 contains special examples, and Section 5 concludes the paper. 



2 THE INTEGRAL REPRESENTATION 

The master equation of our approach is given by the matrix ordinary differential equation 
{iu + H{x)) G{x) = e-^ I , n{x) ^ ( -(^^;HA) ^j 

where x & Mis a real variable; A, u, ^ are real parameters, i = I is unit matrix, and 

G{x) is unknown matrix. Equation (7) appears from the Dyson-Gorkov equation which 
describes two-point correlation functions (Green's functions) of ^Hc — A [26]. Under cer- 
tain physically motivated conventions, Eq.(7) appears from a spatially three-dimensional 
general Dyson-Gorkov equation, and the variable x turns out to be due to a separation 
of the system into a collection of one-dimensional subsystems (spherical coordinates are 
chosen). The parameter A depends, generically, on the angle variables. The parameter 
u (7) corresponds to a thermal Matsubara frequency. Solutions to the Dyson-Gorkov 
equation describe the physics of the model in question [26, 27]. 
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It is convenient to place the Dirac 5-function in the R.H.S. of (7): 

{icu + n{x)) G{x, x') = e'^'^Six -x')l, (8) 

so that (7) appears after integrating (8) over x' from — oo to oo with G{x) = J G{x, x')dx' . 
Let us put (8) in another equivalent form: 



ii-^ + M{x))Gix,x') ^ e'^'^Six-x')!, Mix) = ( . ^ a; + ^A\ 

dx \-x + i/\ -luj J 



where G = cxG, and a = diag{l, —1} is the diagonal matrix. 

Following [26], let us assume that Go{x) is known which respects 



^^G'o(a;) = Go{x)M{x). (10) 



Being so, one is concerned, instead of (9), with solution of 

-^{Go{x)G{x, x')) = -e'^'^Six - x') Go{x) . (11) 
dx I 

Without going into details of [26], let us calculate 

J{x) = JJ e-'''^G{x,x')dx'di 

(.^-integration is also from — oo to oo). Notice that G{x,x') depends on ^ due to (8). 
Integrating (11) one obtains: 

J,,{x) = -(detGo)-' / \Gf{x) r Gl\s)e'^^^-^Us 



Gl^x) r Gf{s)e'^^'-^Ms 

J —oo 



d^, (12) 



where the entries Gq respect (10). Equations (11), (12) tell us that inhomogeneous Eq.(7) 
is solved, in fact, by the method of variation of arbitrary constant [29]. 

In its turn, (10) can be reduced to the couple of second order (Weber's, [3]) equations 
as follows. Let us note a useful matrix relation: 

- '1 = (a- -a-)*-")" "-^^d (13) 

where a = diag{l, —1}, = x ^ d/dx, A = cj + iA, and superscripts and * imply 
transposed and complex conjugated, respectively. Taking Gq in the form 

Goix) = V2(Jj^ j'Ju-\ (14) 

and using 

A* -a-) [X* -a+) ~ [ a-a+ + |A|V ' 
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one obtains unknown /i 2 = fi,2{x) and hi^2 = hi^2{x) from the couple of Weber's equa- 
tions: 

where y = x\^, the tilde implies that /, h are expressed through y, and 1, 2 arc omitted. 
Besides, A is also rescaled in (14'), i.e., A = (c^; + i^)/\/2 until the end of the paper. 

Eventually, the following combinations for /ii 2 and /i^2 are chosen in [26] to express 
G^Q which ensure the integration in (12): 

/ii ± /i - U+{y) T \-^U.{y) , ±i{h T h2) = U+{-y) ^ X-'U_i-y) , 

where U±{y) = W(|A^| ± 1/2, y) are the parabohc cyhnder functions (the notation U{a, y) 
implies -D_(j_i(y)). The asymptotical behaviour olU{a,y) is given by the estimates [3]: 



(2 



,1/2 



W(a, y) ^ exp(-yV4) W(a, -y) ^ ^^^^^ exp(yV4) y^-V^ , 

at y — s> +00, li/l ^ a (the leading terms). 

The theory of distributions allows to calculate (12) straightforwardly. For the given 
choice of Gq, the order of integrations in (12) can be changed at fixed x (Fubini's theorem). 
Integrating over ^ first, we obtain 27r5(s — x). Then, using the properties of 5-function 
[29] we obtain: 



Jiiix) = in 



Gl'{x)Gf{x) + Gl''ix)Gf{x) 



Gh\x)Gf{x) - Gl\x)Gl\x) 

(again in the x-notations). A note about ^^21 will be given in the last section. Other 
entries are not of interest for us. The final answer for J'ii{x) reads: 

Jiiix) = -^(^^Y\xT{\X\' + l)U+{xV2)U+i-xV2) 



+ X*Ti\\\^)U-ixV2)U-i-xV2)) , (15) 
where r(-) is the gamma function [2]. 



3 THE SERIES REPRESENTATION 

Let us consider another approach to (7), (8) [27]. The operator H{x) (7) is self-adjoint 
on M, and it can be conjugated by the constant matrix u (13): 

A ia~ 



-inu - nu oj - ( ia \ 

V — ^a —A / 



Squared operator 7Yg„ resembles the Hamiltonian of a spinning electron in a constant 
homogeneous magnetic field [28]. It is straightforward to obtain eigenvalues Eq, ±En and 
eigenfunctions ^0, {n > 1) oi Hemix): 
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where s = ±, i?„ = a/A^ + 2n , and ijjnix) are the Hermite functions, i/^nix) — tt"-*-/^ 

Hn{x) [30]. OrthonormaUty and completness of the system 
(n > 1) can be directly verified since ipn{x) (n > 0) are orthogonal and complete in 

Let us take unknown G in the form (14) and pass to the equation 



(the first column of the matrix equation in question) . Now the unknown h and / depend 
on X and x', and we expand ( „ ) in ^0(2^); ^n(^) with the coefficients dependent on x': 



f 

■I ^ s=+- n=l 

We define the coefficients [34] and, using 

G\\{x) = j dx'{h{x,x') + f{x,x')), 

obtain J'u in the series form: 

1 TT ^ ,0. . ^ A A 



+ TTJT— ■ (16) 



+ n + 1 |Ap + n 



In [26, 27] one is concerned with (15) and (16) themselves which correspond to the 
integrated entries of G{x) (7). Let us now equate separately the real and imaginary parts 
of (15) and (16). After the analytical continuation |Ap — > z, we deduce: 

r(^)w(^z-^,a;v^)w(^z-^,-xv^) = v^f;^^, Rez>0. (17) 

The series in R.H.S. of (17) is convergent absolutely and uniformly on M, and it is straight- 
forward to express it in the integral form using 

roo 

(n + z)-^ = / e-'^^'+'Ut , Rez>0. 
Jo 

Using the generating function of Hermite polynomials [2] to sum up after exchanging 
summation and integration, one just obtains Eq.(6): 



r{z)u[z-^-,x^ul^z-^,-x^ 



1 r°° f ,1 , x"^ t\ dt _ , , 

= exp U(--2;) tanh - —= , Rez>0,xeR (18.1) 

{x is rescaled, x\/2 — > x). R.H.S. of (18.1) can be transformed to the integral which is 
known [6, 31]: 

il-sY-^ ( x^ \ ds ^fl \^ fl 1 x^' 



t^[\-sY-' \ ds n/1 \^ /I I x^\ 

io (1 + s)^ ^\ 2 ) yf^ V2' ; V2 2' ' 2 y ' 

(18.2) 

where 0i(q;, /3, 7; m, v) is the confluent hypergeometric series in two variables [1, 32]. 
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4 THE SPECIAL CASES 

Two ways of solving (7) give us, naturally, two representations for the same fundamental 
solution G{x). Generally, it is a powerful technique - to equate equivalent representations 
of a given Green function [33]. We shall consider three particular examples in order to 
argue directly the interesting reductions of (18) to the known facts. Equation (18.1) 
depends on two parameters, z eCT and x & M, and their admissible values constitute the 
domain V = {Kez > 0} x M. So, we shall consider the following three subsets in V: the 
line, z = 1/2; open half-plane, x = 0; asymptotic regions, z 1, \x\ ^ is bounded, and 
\z\ is limited, \x\ ^ 1. 

Example 1: z = & M. 

Using [2], we proceed: 



U{0,x)^D_.{x) = [-) K. 



X \ ^/'^ ^, ( x^ 



In their turn. 



2^ 



Then, 



and 



x"'M(0, .)M(0, -X) = ^ (^) (/. (t^ + /_ . y (19) 



On the other hand, the L.H.S. of (18.2) can be written at ^ = | as 



/' 

JO 



'J? \ ds 



Accordingly to [31] (Chapter 2), the integral (20) is just equal to (19) (see APPENDIX). 

Example 2: Re2; > 0,a: = 0. 
At X = the L.H.S. of (18.1) is [2]: 

r(.)w^(.-i,o) ^ r{z)DlM - 
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R.H.S. of (18.1) takes at a; = the form of the definition of the beta function [2]: 

i|../-(l-.)-V^d.^iB(ii). Re(|)>0. 
The coincidence is clear. 

Example 3: z ^ 1, \x\ 7^ is bounded, or \z\ is bounded, \x\ ^ 1. 

In leading order, we estimate the triple product r{z)U{z — |, x)U{z — |, —x) in the L.H.S. 
of (18.1) as [3]: 

r(.) « (. - 1, .) « (. - i, -.) . ( j,^) , (21) 

where z > 1/2, and + 4^; — 2 is assumed positive and large. By steepest descent we 
estimate the R.H.S. of (18.1) at 2; » 1 (|a;| 7^ is bounded) as / {2z))^^'^ . Analogously, 
the integral in the L.H.S. of (18.2) can be estimated at \x\ » 1 (|2;| is bounded) as 
(27r)^/^/|a;|. In the both cases we get agreement with (21). However, when z and \x\ are 
large, (18.1) results in (21) as well. 

5 CONCLUSION 

To conclude, the following remarks are in order. The relation 

D,{z) = 2^+2^1^., 1^1 

^ ' y/z 2 + 4'±4 y2 ) 

is known which relates Di,{z) to the Whittaker function. Thus, our product Di,{z)D^{—z) 
can be expressed by means of the integral representation for iyK^^(z)VF^,^(z') [2] (though 
containing the hypergeometric function 2F1 under integration). However, some extra work 
is still required to relate (6) to the corresponding formula. Consideration of the entry J'21 
does not give really different formula: we obtain a representation which can be deduced 
if one applies x — d/dx to (6) (compare with (3)). The representation (6) looks simpler 
in comparison with analogous formulas mentioned above, and hopefully can be useful in 
applications. 
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APPENDIX 

More precisely, [31] provides the formula for the integral 

r ^""^ -px > 

Jo \/o? — x^ 



8 



which can be obtained from 



e-P'^dx (7rj9)i/2 

— Ki 



ap 

q2 _ j.2\ 2 * V 2 



X (^i (f ) + ^-i (f)) ' a>0,Rep>0 (a) 
by multiple differentiation p)". Besides, there exists another integral formula [31]: 



/ 

J a 



— K 



^Jx{x^ — a^) 



2 3 V 2 



^sbxd^ ^ fnb\ kJ^-^] I J'-^] , b>0,Rez>0. 



X (f ) - (f )) , a>0,Rep>0. (b) 

The couple of Eqs.(a) and (b) can be related to the following couple of the integral 
transforms [35]: 

Jo ^x{x^ + z^) V 2 y H 2 ; 4 2 ; ' 

Notes added after release of the paper 
1. 

It is appropriate to point out that the representation 

exp(-x') (hI{x) + Gl{x)) = 2^+V-^r(A+l) 

X j exp (-(2A + l)t + x^tanhij (coshi sinhij-^/^di (A^l) 

has been found in [36] . Here, Hx and G\ are the Hermite functions which are related to 
the confluent hypergeometric function $(a, 6; z), 

^[a,b,z) = — — 2^ 



r(a) tonb + n) n!' 
as follows [37]: 

^ Arr / N r(l/2) A 1 2n r(-l/2) A I 3 2^ 

^ / X r(l/2) ttA^, a 1 2x r(-l/2) ttA^, A 13 2x 

2 ''^l^' = - r(l/2 - A/2) T 2 ■ 2- - ' + fRM " T 2 + 2- 2' ^ ' ' 

Equation (Nl) appears in [36] as a confluent limit of an integral representation for product 
of two Gegenbauer functions. The representation (Nl) turns out to be useful [36] in the 
problem of monotonicity of differences of zeros of the Hermite functions [38, 39]. 
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Using connection [37] between $(a, b; z) and the parabolic cylinder functions, one finds: 

exp(-xV2) Hx{x) = 2^I^Dx{xV2) , 

(iV2) 

2A/2 

ex.p(-x'^/2)Gx(x) = (cosnXDx(xV2) - Dx(-xV2)) . 

smnX ^ ' 

Eventually, using (N2) and the linear relations [2] between the parabolic cyhnder func- 
tions, it is straightforward to check that (Nl) and (6) (and thus (18)) are equivalent. 



Equation analogous to (17) can also be written for the product of two parabolic cylinder 
functions at coinciding arguments: 

T{z)U'[z-\,x^ = ^p^tll^M^ Rez>0, (7V3) 

where summation in R.H.S. goes over simple poles of r(^). In the same way as in Sec. 3, 
we obtain from (N3): 

D'_,{x) ^ 2-'/'r-\z) /^°°exp((i-z)i-^coth0^£=, Re^>0, (m) 

or 

However, the representation (N5) is known [31] just for the square of the parabolic cylinder 
function but in another form: i.e., for the integration variable obtained by shift s — 1 — > s. 
The integral in (18.2) and Eq.(N5) imply the following combinations: 



POD 

D^,(x)(D_,(x)+D_,(-x)) = r-\z) / 

Jo 



s — I X \ ds 



(JV6) 

- = r-(.-) f 1^ j^exp (4.) I . 

Besides, (18.1) and (N4) can be re- written together in the following form: 

D^.(.)D^.(±.) = i^|(=i|) f exp ((i - .)t ^ .He' ^^r')^ - ^ 

The Author is grateful to Prof. M. E. Muldoon for attracting his attention to [36] 
and for interesting correspondence which inspired to write down Eqs.(N3)-(N7) and to 
present the given notes. 
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